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Abstract 

I This is a brief review for particle physicists on cosmological impact of neutri- 

fT") ' nos and on restrictions on neutrino properties from cosmology. The paper includes 

discussion of upper bounds on neutrino mass and possible ways to relax them, 
1— 1, methods to observe the cosmic neutrino background, bounds on the cosmological 

'"q! lepton asymmetry which are strongly improved by neutrino oscillations, cosmolog- 

1 ■ ical effects of breaking of spin-statistics theorem for neutrinos, bounds on mixing 

parameters of active and possible sterile neutrinos with the account of active neu- 
r~| \ trino oscillations, bounds on right-handed currents and neutrino magnetic moments, 

and some more. 

(N : 

> : 

£j : 1 Introduction 

00 ■ 

Neutrino is the weakest of all known elementary particles, but despite that, or maybe 
CO because of that, cosmological impact of neutrinos is significant. Neutrinos are important 

for cosmology and, vice versa, astronomy allows to measure neutrino properties with 
O ■ precision which is in many cases higher than the precision of direct experiments. 

J> ! Neutrinos are the second most abundant particle in the universe (after photons of 

I CMBR), their total number density, with antineutrino included, is 

S i K + n p )^ 340/cm 3 , (1) 

almost equally shared between all three active neutrino species u e , u T , and v T . For com- 
parison the average cosmological number density of the normal baryonic matter is 

n B « 2.5 x 10- 7 /cm 3 . (2) 

The average momentum of cosmic neutrinos is very low, about 6.1° K ~ 5.3 ■ 10 -4 eV. The 
cosmic neutrino background (Cz/B) is not observable directly at the present time, because 
of the very weak interactions of these low energy neutrinos. Still their cosmological impact 
is profound. 
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In this talk I will consider cosmological effects produced by neutrinos and cosmological 
bounds on neutrino properties. More detailed discussion can be found in the review 
paper [T], but presentation here contains some recent development. 

The content of the talk is the following: 

1. Cosmological bounds on neutrino masses, m„. 

2. Possibility of direct registration of C^B. 

3. Restriction on the number of neutrino species from BBN and CMBR. 

4. Neutrino statistics and cosmology. 

5. Neutrino oscillations, BBN, and cosmological lepton asymmetry. 

6. Right handed currents and magnetic moment of neutrinos. 



2 Thermal history of neutrinos and Gerstein-Zeldovich 
bound 

In the early universe at T > 1 MeV neutrinos were in thermal equilibrium with electron- 
positron pairs and photons. Correspondingly their number density was 

n u = (3/8)n 7 (3) 

for each left-handed neutrino flavor, u e , i/^, u T , under assumption of vanishing asymmetry 
between neutrinos and antineutrinos, i.e. 

n u = n p . (4) 

The condition for thermal equilibrium is found by comparing the neutrino reaction rate, 
r = an ~ G 2 F T 5 , with the cosmological expansion rate, H = a /a ~ T 2 /mpi. If T > H, 
thermal equilibrium is established. This took place at T > 1 MeV. 

More accurate treatment based on solution of the kinetic equation governing neutrino 
distribution pQ shows that v e decoupled from the electron-positron pairs at T m 2 MeV, 
while v^ T decoupled at T m 3 MeV. Interactions between neutrinos maintain their kinetic 
equilibrium down to lower temperatures, T Ve = 1.3 MeV and T„ VjVr = 1.5 MeV for v e 
and v^ T respectively. Below these temperatures neutrinos practically decoupled from the 
plasma and free streamed with the speed of light till they became nonrelativistic. 

At smaller temperatures, T < m e , photons were heated by e + e~-annihilation, while v 
were not, because they had already decoupled from the electromagnetic component of the 
primordial plasma. The increased number density of photons had led to the decrease of 
the neutrino-to-photon ratio: 

n v + n p = (3/ll)ra 7 = 112/cm 3 (5) 

The coefficient 3/11 is obtained from the entropy conservation and the numerical value of 
the neutrino+antineutrino number density is presented for the present day universe and 
obtained from the measured number density of photons in cosmic microwave background 
radiation (CMBR), n 7 = 410.5 ± 0.5/cm 3 . 
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As a result of the photon heating the temperature of neutrinos dropped down with 
respect to the photon temperature: 



T v = 0.714T 7 = 1.945AT" = 1.68 • l(T 4 eV, (6) 

where the second equality is written for the present day neutrino temperature calculated 
from the known value of T gamma = 2.725/^°. One should remember, however, that for 
massive neutrinos the distribution function has the form: 

/„ = [exp(p/T) + l]" 1 ■ (7) 

It is not equilibrium distribution because in the equilibrium one there should be the 
neutrino energy, E, instead of momentum, p, as in eq. (CO). So, strictly speaking, parameter 
T here is not temperature. The difference, however, is essential only at low temperatures, 
T ~ m u . 

Knowing the number density of neutrinos in the present day universe we can calculate 
their energy density and from the condition that the latter does not exceed the measured 
energy density of matter we obtain the cosmological bound on neutrino mass: 

E a m„ a < 94fi TO /i 2 eV (8) 

where h = H /100km/ s/Mpc is the dimensionless Hubble parameter and Q m is the fraction 
of the energy of matter relative the critical energy density, p c = 3H 2 m 2 Pl /87T. This bound 
was derived by Gerstein ans Zeldovich in 1967 [2J. 

Using the observational data: Q m < 0.25 and h = 0.7, we find: 

E a m Va < 11.5 eV (9) 

For almost equal masses of neutrinos (as we know from neutrino oscillations): 

m„ a <3.9eV (10) 

This bound can be further improved by factor ~3 because too large energy density of 
neutrinos would inhibit large scale structure (LSS) formation at relatively small scales 
and at hight redshifts, z > 1. In simple words, neutrinos are fast and reluctant to 
form gravitationally bound systems. To eliminate such an undesirable property the mass 
density of hot dark matter, created by cosmic neutrinos should not exceed ~30% of the 
total matter density. Correspondingly: 

m Ua < 1.2 eV (11) 

This is a robust and quite strong bound, stronger than that obtained in the direct exper- 
iments. 

The best bounds on neutrino mass, obtained from tritium decay experiments (Troitsk 
and Meinz) is [I]: 

m Ue < (2 - 3) eV (12) 
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Though this bound is usually presented as a bound on the mass of electronic neutrino, it is 
inaccurate because, as we know from the measured neutrino oscillations, mass eigenstates 
are strongly different from the flavor eigenstates, i.e. from v e ^ T . Neutrino oscillations 
allow to measure only two mass differences of three mass eigenstates, but not the absolute 
value of the mass. Neutrino oscillation data are best fitted with 

5m 2 solar = (5.4 - 9.5) ■ 1(T 5 eV 2 

5m 2 tm = (1.2-4.8)-10- 3 eV 2 (13) 

Presently telescopes allow to weight neutrinos more accurately than direct experiments. 
Astronomy can be sensitive to 

m u ~ (a few)x0.1eV, (14) 

based on combined data on LSS and CMBR (see sec. EJ). 

Historical remark: the paper by Gerstein and Zeldovich "Rest mass of muonic neutrino 
and cosmology" was published in 1966 [2]. Six years later a similar paper by Cowsik and 
McClelland j3] "An upper limit on the neutrino rest mass" was published. In many sub- 
sequent works the cosmological bound on neutrino mass is quoted as "Cowsik-McClelland 
bound". This is not just, however, firstly, because the GZ paper was much earlier and, 
secondly, in the paper by Cowsik and McClelland the effect of photon heating by e + e~- 
annihilation was disregarded and both helicity states of neutrinos, left-handed and right- 
handed, were assumed to be equally populated. This incorrectly gives rise to 7 times 
stronger bound. 

We know that only left-handed neutrinos participate in weak interactions. If neu- 
trino mass is nonzero the other, right-handed state, must exist. However, right-handed 
neutrinos could be in equilibrium only at very high temperatures and even if the were 
abundantly created at some early cosmological stage, they would be strongly diluted 
by entropy released in massive particle annihilation and so their number density would 
normally be negligible at T ~ MeV. 

2.1 Is it possible to relax the GZ bound? 

Let us critically discuss essential assumptions used in the derivation of the GZ bound and 
check if they can be modified in such a way that the present day number density of CiaB 
would be noticeably smaller than the standard one ([5]) and thus would allow for a larger 
neutrino mass. 

The obtained upper bound on m v is based on the following: 

1. Thermal equilibrium in the early universe. It is surely true in the standard cosmology. 
Neutrinos would not be produced in equilibrium amount if the primeval plasma temper- 
ature were never larger than a few MeV [5]. It could be realized in inflationary models 
with anomalously low (re) heating temperature. However, a smaller number density of 
neutrinos during big bang nucleosynthesis (BBN) would distort successful predictions for 
light element abundances. The problems with baryogenesis would be also serious. 

2. Non-vanishing lepton asymmetry. We assumed that number density of neutrinos is 
equal to that of antineutrinos. If the primeval plasma has a non-zero lepton asymmetry 
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this equality would be broken. However, in thermal equilibrium the total number density 
n u +n p would be larger in asymmetric case and the upper bound on m v would be stronger. 
Moreover, as we see below, the cosmological lepton asymmetry is strongly bounded from 
above and cannot noticeably change the neutrino number density. 

3. Conservation of the ratio n v jn^ from 1 MeV to the present day. If there are some extra 
sources of heating of the photon plasma below MeV, the ratio n„/ri 7 would be smaller 
than (jSJ). This can be achieved e.g. by electromagnetic decays of new light long-lived 
particles. Care should be taken of BBN and the frequency and angular spectra of CMBR. 
In particular, the ratio %/n 7 measured at BBN and CMBR cannot differ more than by 
factor 2. It seems to be difficult and one has to invoke some more new physics, but maybe 
this is not excluded. 

4. Neutrino stability at the cosmological time scale. It is assumed that neutrinos created 
in the plasma at T ~ 1 MeV survived up to the present time, i.e. r u >t u ~ 10 10 years. 
Bearing in mind a small neutrino mass (1121) . it is difficult to avoid the conclusion about 
a long life time of v. To facilitate neutrino decay new interactions are necessary. Even 
if such exotics exists, it does not lead to relaxation of GZ limit if neutrino decays into a 
very light or massless particle (e.g. photon) plus another a little lighter neutrino, because 
in this case the total number of neutrinos is conserved. As we see from eq. ([TBI the mass 
difference is much smaller than eV and thus the upper bound on the neutrino mass re- 
mains the same. However, if the decay goes into a new sterile neutrino, the cosmological 
bound on the masses of normal active neutrinos could be weaker. 

5. An absence of new strong interactions between neutrinos at low energies. A new inter- 
action may enhance //^-annihilation into new very light particles, so the neutrino number 
density today would drop down. To this end a very light or massless boson is needed. 
Such bosons are practically excluded by stellar cooling, see e.g. [6]. 

6. Absence of right-handed neutrinos. If both helicity states of neutrinos would be popu- 
lated the bound would be stronger. 

7. Presented below stronger bounds, based on angular spectrum of CMBR and analy- 
sis of LSS, depend upon the values of the cosmological parameters, on the spectrum of 
primordial perturbations, and on the type of dark matter. A nonstandard spectrum of 
density perturbations at small scales may change these strong bounds on m v but would 
not influence the classical GZ bound. 

2.2 Non-equilibrium cosmic neutrinos 

Usually the cosmological expansion does not destroy equilibrium distribution of massless 
particles, even after their interaction is effectively switched-off. We observe that by perfect 
black body spectrum of CMBR. However, this is not exactly true for neutrinos. As 
we mentioned above, neutrinos decoupled from cosmic plasma at T = 2 — 3 MeV. At 
that moment the primeval plasma consists of four weakly coupled parts: electromagnetic 
(e + ,e~,7) and three separate neutrino parts. At T ~ MeV the temperature of photons 
became larger than the temperature of neutrinos due to the heating of photons by e + e~- 
annihilation. Residual annihilation of hotter e + e~-pairs into neutrinos heats up neutrinos 
and distort their spectrum [7j. For more details and references see review pQ. This leads 
to a larger neutrino number density, which can be described as an increase on the effective 
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equilibrium number of neutrino species: 



AN, 



V 



0.03 + 0.01 



(15) 



Here the second term, 0.01, comes from the plasma corrections which diminish n 7 with 
respect to the ideal gas approximation [8]. Because of these corrections the number in 
eq. ([T]) is not 336/cm 3 , i.e. three times of eq. but slighty larger. 

This correction has a negligible impact on primordial 4 i7e abundance, because there 
are two effects of the opposite sign. An increase of the effective number of neutrino species, 
AiV > 0, leads to an increase of the mass fraction of the produced helium. On the other 
hand, the account of the spectrum distortion of v e and their larger number density has 
an opposite effect on A He and the net effect is quite small. 

However a larger energy density of the relic neutrinos has an impact on the form of 
the angular spectrum of CMBR which may be noticeable in the future CMBR missions, 
in particular, in Planck. The observation of this effect will be a measurement of physical 
processes which took place, when the universe was only 1 second old. 

3 Neutrino and large scale structure formation 

Investigation of the large scale structure of the universe is an essential ingredient of the 
modern cosmology. Comparison of the theory with observations allows to study the prop- 
erties of dark matter, measure cosmological parameters and the spectrum of primordial 
density perturbations. In particular, an analysis of LSS permits to put more stringent 
than (FlOi) bounds on neutrino mass. Bound (1111) is already based on the consideration of 
LSS but a more detailed study leads to stronger bounds. 

3.1 Basics of LSS formation theory 

Theory of LSS formation is based on the following input: 

1. Spectrum of primordial density fluctuations. The Fourier transform of the power 
spectrum, i.e. of {b~p) 2 is usually parametrized in the simple form: 



Observational data indicate that n ~ 1, i.e. the spectrum is close to the so called fiat, or 
Harrison-Zeldovich, spectrum. For such spectrum metric fluctuations, which are dimen- 
sionless, do not contain any dimensional parameter, 



This type of spectrum is predicted by inflation with n slightly different from unity. 

At large scales, from Gpc down to about 10 Mpc, the spectrum is measured by the 
angular fluctuations of CMBR. At the other end, from Mpc up to tens Mpc the spectrum 
is measured by the observed LSS. A good test of the theory is that these two kinds of 




(16) 




(17) 
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independent measurements give coinciding results in the common region of wave lengths 
around 10 Mpc. 

2. Type of the density perturbations. It is assumed that the density perturbations are 
the so called adiabatic ones, i.e. they can be visualized as time-shifted in one space point 
with respect to another, p = p(t(x)), as predicted by inflation. In other words, density 
perturbations are the same in all forms of matter: baryons, dark matter, photons, etc. 
This type of perturbations is confirmed by CMBR at large scales, above 10 Mpc. It is not 
excluded that at smaller scales the so called isocurvature perturbations may be significant, 
but usually they are neglected. 

Comment: density perturbations at the early stage can be decomposed into two in- 
dependent modes: adiabatic and isocurvature. We have already defined adiabatic. The 
isocurvature ones can be understood as perturbations in the chemical content of the pri- 
mordial matter and initially they have vanishing energy density perturbation, 5p = 0. 
Later, when the equation of state becomes different due to different matter content the 
isocurvature perturbations give rise to the normal density perturbations, 5p ^ 0. 

3. Properties of dark matter. According to observations the total energy density of the 
universe is close to the critical one, 

p^ Pc = „ io-2 9 g/cm 3 (18) 

In terms of the dimensionless cosmo logical parameter Qj = Pj/p c we have: 

Q tot = 1 ± 0.05, tt B ^0.05, n DM ^0.25, Vt DE «0.7. (19) 

where B, DM, and DE are respectively baryons, dark matter, and dark energy. 

Usually it is assumed that the cosmo logical dark matter consists of practically non- 
interacting cold dark matter (CDM) particles. Dark energy remains mysterious but ac- 
cording to observations its equation of state is quite close to the equation of state of 
vacuum energy, p = —p, and hence dark energy is either vacuum energy (i.e. Lambda- 
term) or something similar to it. So the standard model is called CDMA model. Still it 
is not excluded that dark matter may be self-interacting as e.g. mirror matter and dark 
energy may have equation of state which varies with time, p = w(t)p. 

3. Analytical calculations at linear regime, when Sp/p <^ 1. The calculations are done in 
the frameworks of classical physics: general relativity and hydrodynamics. The results 
are not distorted by subsequent evolution at large scales, larger than 10 Mpc, accessible 
to CMBR. 

4. Numerical simulations at nonlinear regime, when 5p/p > 1. The simulations are neces- 
sary at smaller scales, < 10 Mpc. Simulations are usually done with dark matter particles 
with masses of about 1O 6 M , where M Q = 2 • 10 33 g is the solar mass. Some essential 
physics in the simulations may be omitted especially at smaller scales but the overall 
picture seems to be fine. 

If any of the above is not true, the strong bound on m v may be relaxed. On the other 
hand, all above is confirmed by the data. One cannot say that all the assumptions are 
strictly proven, so some freedom remains, but the agreement of this simple picture with 
the observations is quite impressive. 
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3.2 Neutrino role in LSS formation 



Neutrinos stopped to interact with the primeval plasma at T ^> m v i.e. they were 
relativistic at decoupling. After that they freely propagate in the universe with the red- 
shifted momentum. Their free streaming length before they become nonrelativistic is 
huge: 

l FS « 250 Mpc (eV/m v ) (20) 
The mass inside the free-streaming volume is 

M FS = « °- l! ^T « 10 17 M Q (eV/m,) 2 (21) 

In neutrino dominated universe the structures below this scale should be erased. If neu- 
trinos are sub-dominant they are still able to inhibit the structure formation at small 
scales. 

In the universe dominated by neutrinos very large structures with the size of the 
order of Ifs should be formed first and after that smaller structures could be formed 
by fragmentation of large scale structures to smaller ones. In other words smaller size 
structures are younger. The same effect but less pronounced would exist in mixed hot-cold 
dark matter scenario. Thus the larger is the neutrino mass (but still small) the longer 
would be the delay in the formation of small structures. In particular, the number of 
Ly-a clouds at high red-shifts, z > 1, is very sensitive to light massive neutrinos. The 
data analysis permits to restrict the neutrino mass by: 

£m„<0.7eV, (22) 

i.e. the individual neutrino mass is bounded by m v < 0.23 eV. 

On the other hand, formation of Lyman-a clouds is sensitive to gasodynamics and 
possible shock waves which are not yet included into numerical simulation at the non- 
linear stage of perturbation evolution and are poorly understood. Another unresolved 
problem is that of biasing, i.e. of the shift between the distributions of the visible and 
dark matter. The biasing is observationally determined at galactic scales or larger and is 
extrapolated to smaller scales but it is unclear how much it changes with the size of the 
object. 

An ad hoc modification of perturbation spectrum at intermediate scales could mimic 
presence of massive neutrinos and a less restrictive mass bound from LSS would be ob- 
tained. 

On the other hand, a good agreement between perturbations at the scales accessible 
to the measurements of angular spectrum of CMBR, especially by WMAP, and smaller 
scales, measured in the analysis of LSS by 2dFGRS and SDSS, in the common region, 
accessible to both CMBR and LSS, demonstrates that the picture is self-consistent. 

For a review of different methods of determination/restriction of the neutrino mass 
from the astronomical data see ref. |10j . 
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4 Neutrino mass and CMBR 



The spectrum of the angular fluctuations of CMBR provides us with accurate information 
on cosmological parameters and, in particular, on neutrino mass and on the number of 
neutrino species. Theoretically calculated spectrum of CMBR angular fluctuations for the 
accepted values of the cosmological parameters is presented in fig. 1. The shape of the 
spectrum depends upon neutrino mass and its analysis allows to put a stringent bound on 
the sum of the masses of all three neutrino species. In fig. 1, taken from ref. [9] five curves 
corresponding to different neutrino masses m u = 0.5, 0.7, 1.0, 2.0 eV are presented and 
compared with the observational data from WMAP3. 

Massive neutrinos lead to the following two effects in the angular fluctuations of 
CMBR: 

1. Shift of peaks to the left with rising mass. The larger m u , the earlier is the non- 
relativistic stage. Thus the distance to the last scattering surface is shorter and peaks 
move to smaller I. However, this effect can be compensated by a shift in H. 

2. Decrease of the (1st) peak hight. Neutrinos with m„ > 0.6 eV become nonrelativistic 
before recombination, so the matter radiation equality takes place earlier and the en- 
hancement of the first peak by the integrated Sachs- Wolfe (ISW) effect becomes weaker. 
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Figure 1. Effects of massive neutrinos on the spectrum of angular fluctuations of CMBR. 
The curves from top to bottom in the first maximum correspond to vn v = 0, 0.5, 0.7, 1.0, 
and 2.0 eV. 

The degeneracy of the effects induced by the neutrino mass and the magnitude of the 
Hubble parameter is illustrated in fig. 2 [5]. On the same figure the measured values of 
H by three different groups are presented. 
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The analysis of WMAP data allows to restrict the neutrino mass from above by [9]: 



m v < 0.63 eV 



(23) 



at 95% confindence level. The derivation does not use data on LSS which are subject to 
ambiguities due to poorly understood biasing and nonlinearity in structure evolution. 
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Figure 2. The constraints on h for several fixed values of neutrino mass. The first three 
bars at the left hand side are the constraints from distance ladder measurements. 



5 Massive conclusion 

Here we summarize 

1. All bounds quoted above are based on the calculations of the present day number 
density of neutrinos initiated by GZ. 

2. From CMBR only it follows: m v < 0.63 eV. 

3. More restrictive bounds obtained in several papers under different assumptions vary 
in the interval: Em„ < 1 — 0.42 eV. These results are based on the analysis of LSS at 
relatively small scales plus CMBR to shift degeneracy of parameters. 

4. With almost equal neutrino masses one can conclude: m v < 0.63 — 0.14 eV, lower value 
is probably too strong. 

4. The planned detector KATRIN, expects to reach the accuracy m u < 0.35 eV at 5a. 
Do we need this expensive experiment? 

5. Possible ways to violate the upper limits on m v : 
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a) From particle physics: new particles or stronger than standard ^-interactions; 

b) From cosmology: unusual spectrum of (isocurvature) perturbations at small scales or 
low (re)heating temperature. 

So to conclude, we need KATRIN. The issue of neutrino mass is too important to rely 
only on astronomical measurements. However, it is quite probable that the absolute value 
of the neutrino masses is close to their mass differences. In this unfortunate case KATRIN 
will not help but there may be a chance to weight neutrinos by the future Planck mission. 

6 Detection of cosmic neutrino background (Cz/B) 

Though the number density of cosmic neutrinos is very high, it is not yet directly reg- 
istered. There are some indirect ways to "observe" Cz/B, for example by its impact on 
LSS and CMBR in the case of sufficiently large m u , as is discussed above. Among other 
possibilities there are: 

1. Measurement of the number of neutrino families through BBN. Larger number of equi- 
librium neutrino species shifts n/p-freezing to higher T and shortens the time of nuclear 
formation. One neutrino family makes 4 He 5% larger. With ng/ri-y known from CMBR: 

N v = 3± 0.3 (24) 

(for more detail see below sec. [7]). 

2. The angular spectrum of CMBR. A change of the number and energy density of the relic 
neutrinos would shift the moment of matter-radiation equality, i.e. the moment when the 
universe stopped to be dominated by relativistic matter and the cosmological expansion 
regime changed into the matter dominated one. This would shift the positions of peaks of 
CMBR spectrum. The effect is rather weak and the precision is lower than that of BBN. 
For massless or very light neutrinos the number of effective neutrino families is [UJ 

N v = 3 ± 1 (25) 

It is clearly a long way to the necessary level of accuracy to observe the effect discussed 
in sec. E21 eq. (fT5l). 

3. Direct detection of CvB: 

a) The evident process of registration of the relic neutrinos by i/A r -scattering seems to be 
out of question because the corresponding cross-section is extremely small: 

o vN ~G 2 F El~\^cm 2 (J^f (26) 

Since the average energy of massless cosmic v is (E) m 3.15T m 5.3 • 10~ 4 eV we obtain 
a V N ~ 10 -62 cm 2 . Even for massive neutrinos the situation is not much better. 

b) Z-burst effect [12] opens a more realistic chance to register Cz^B. To this end a flux 
of ultrahigh energy cosmic ray neutrinos is necessary. The Z-resonance scattering of the 
high energy neutrinos on the background ones has much larger cross-section: 

<j z ~ a/m% « 1(T 33 - 1(T 34 cm 2 (27) 
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However, to excite the Z-resonance on massless v with E ~ 10 3 eV very high energy of 
cosmic ray neutrinos is necessary: 

E v , hi = 4 • 10 21 eV (eV/E U]low ) (28) 

Even if neutrinos with E ~ 10 24 — 10 25 eV exist in the cosmic rays, their flux most prob- 
ably is negligibly small and the chance to register the Z-burst is practically zero, but it 
may be feasible for massive, not too light v. 

c) The cross-section of elastic scattering of low momentum Dirac neutrinos can be strongly 
enhanced by coherent effects. Neutrinos efectively interact with all nuclei inside the vol- 
ume A 3 , where A = l/m u v is the inverse neutrino momentum. Correspondingly the 
cross-section would be larger than that for a single nucleus by the factor iV 2 , where N is 
the number of nuclei inside A 3 . According to the estimates of ref. [13] the acceleration of 
a small piece of matter with the size I ~ A could reach the value 10 -22 cm/s 2 . This result 
was obtained for m v = 30 eV and the galactic density of neutrinos n^ ga ^ = 10 7 /cm 3 . In 
other words, it was assumed that neutrinos make all dark matter in the galaxy. Now 
we know that neutrinos are much lighter and can make at most 15% of all dark matter. 
Correspondingly the effect should be 6 times weaker or maybe even more than that if light 
neutrinos are not accumulated in the galaxy in the same fraction as cold dark matter. 

d) An interesting possibility is open by the inverse beta decay, induced by the relic neu- 
trinos: 



v 



A^A' + e, (29) 



This process should take place if nucleus A /^-decays into nucleus A'. The reaction has zero 
threshold but if m„ = the spectrum of the produced electrons in the inverse beta decay 
very weakly deviates from the spectrum of electrons in the beta-decay. As was suggested 
in ref. [H], the spectrum distortion may be observed if the chemical potential of CzaB 
is large, fi u > T v . However, BBN plus large mixing angle solution to the solar neutrino 
anomaly strongly restricts chemical potential of all neutrino flavors [15] (see sec. [7]): 

fi u < 0.077; (30) 

The /3-capture of massive neutrinos P3] looks more promising for registration of CzaB, if 
m v is not too small. The allowed kinematical region for the electron energy in /3-decay is 
m e < E e < Wo, while in the //-capture it is E e > W + 2m u . According to the calculations 
made in ref. [TB], the ratio of the reaction rates for the detector resolution A is equal to: 

A/3 (A) _ 2 fA\ 3 A + 2m,\ 3/2 >> 1 



K 9C(3) \TJ 

where T v « 1.7- 10~ 4 eV. 

The above estimate ( 131]) is valid for the conventional Cz/B with the distribution func- 
tion given by eq. ([7|) but if by some reason there are additional contributions into the 
number density of the cosmic neutrinos the result may be quite different. In this connec- 
tion is tempting to ask if the events observed in experiment [T7] out of the kinematically 
allowed region of the tritium beta-decay be explained by some anomalous contribution 
into f v l 
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According to the calculations of ref. [16] the neutrino capture rates in events/year, 
for 100 g of 3 H, for the Fermi-Dirac distribution (J7|) over neutrino momentum but not 
energy, with T„ = 1.7 • 10~ 4 eV (FD), for a Navarro- Frenk- White profile (NFW), and for 
the present day mass distribution of the Milky Way (MW) are respectively: 7.5 (FD), 
23 (NFW), and 33 (MW) if m v = 0.3 eV and 7.5 (FD), 10 (NFW), and 12 (MW) if 
m v = 0.15 eV. 



7 Neutrino and big bang nucleosynthesis (BBN) 

Big bang nucleosynthesis took place when the universe was quite young, t = 1 — 300 sec, 
and not too hot, T = 1 — 0.07 MeV. At this period primordial light elements, 2 H, 3 He, 
4 ife, and 7 Li, were created. Physics in this energy range is pretty well known and the 
results of calculation are in a good agreement with the observations. Deviations from the 
usual physics would destroy the agreement between theory and the observational data, so 
BBN serves as an efficient "cleaner" of distortions of the standard physics. 

The building blocks for formation of light elements were prepared by the reaction of 
neutron-proton transformations: 

n + e + <-> p + u e 

n + v e <-> p + e~ (32) 

These reactions were frozen at Tf ~ 0.65 MeV, and this determines the starting value 
of n/p-ratio, n/p « exp(— Am np /Tf) m 0.135. When T drops down to 60-70 keV (it 
corresponds roughly to 200 sec) all neutrons, which survived the decay, quickly form 4 ife 
(about 25% by mass) and a little 2 H (3 x 10~ 5 by number), 3 He (similar to 2 H) and 7 Li 
(io- 9 - 10~ 10 ). 

It is already clear from reactions ( 1321) that neutrinos play an important role in BBN. 
There are several neutrino effects which influence the primordial abundances of light 
elements: 

1. Impact of neutrinos on the cosmological cooling rate. The energy density of the primeval 
plasma is: 

< &m 2 Pl it 2 . 
P= 32^ 2 = 30 g * T (33) 

where 



g* = 10.75 + 1.75AN U (34) 

is the number of the degrees of freedom of particles which were present in the plasma. 
The first term, 10.75, comes from photons, e ± -pairs, and three families of left-handed 
neutrinos. All the rest (which is absent in the standard model) is parametrized as AN U . 
This AN U includes any form of energy present at BBN. AiV = 1 corresponds to energy 
density of equilibrium neutrinos and antineutrinos with vanishing chemical potential and 
one polarization state. In particular, if neutrinos are degenerate, i.e. n u ^ n p and thus 
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their chemical potential is non-zero their energy density would be larger than the usual 
equilibrium one by: 



AN 



V 



15 U 




) 



2 



(35) 



where £ = ji/T. 

Positive AN U leads: 

1) to an earlier n/p-freezing and higher frozen n/p-ratio; 

2) to a faster cooling down to the BBN temperature, Tbbn ~ 65 keV, and to higher n/p 
as well, because less neutrons would have chance to decay. The net result is an increase 
of A He by about 5%. 

Negative AN leads to sign-opposite results. From the observed abundances of 4 He and 
known from CMBR ratio of baryons to photons, n^/n^ = 6 • 10~ 10 one can conclude that 
AN V < 0.3, though, strangely, AN U < or N v < 3 seems most favored. A recent analysis 
of determination of ris/n^ and AiV from BBN and CMBR can be found in ref. [T8] . 
2. Impact of degeneracy of electronic neutrinos. If chemical potential of v e is non-zero, 
its effect on BBN is much stronger than the effect of comparable chemical potentials of 

and v T . A non-zero £ e shifts equilibrium value of n/p-ratio as: 



Hence the bounds on the chemical potentials of v e and i/^ T are very much different: 



They are valid if a compensation between the effects induced by £ M T and £ e is allowed. In 
the absence of such compensation the bounds are somewhat stronger. We will see below, 
however, that the bounds on chemical potentials of all flavors are much stronger because 
of efficient neutrino oscillations in the early universe. 

3. Variation of the energy density of electronic neutrinos. If the asymmetry between 
neutrinos and antineutrinos is negligible but the energy density of v e differs from the 
equilibrium one, the effect of that would have the opposite sign with respect to the effect 
induced by non-zero AN V . If n Ve > n eq , the n — p freezing temperature would go down 
and so would the n/p-ratio. In the opposite case, if n Ve < n eq the ratio n/p would rise. 

4. The effect on BBN of the oscillations between active neutrinos in equilibrium is neg- 
ligible because in this situation the oscillations do not change anything in the plasma. 
However, if neutrino asymmetries are nonvanishing, the oscillations, which do not respect 
conservation of the flavor lepton numbers, lead to transitions between electronic, muonic, 
and tauonic asymmetries. We know that the large mixing angle solution to Solar neutrino 
deficit is realized. Hence initially large £ M or £ T would be effectively transformed into £ e , 
prior to BBN and all chemical potential would be equalized. Since £ e is strongly restricted, 
similar bounds can be imposed on all i/-flavors [I~5] : 




(36) 



|^r| < 2.5 
|£ e |<0.1. 



(37) 
(38) 



16*1 < 0.07 



(39) 
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Evolution of lepton charge asymmetry for LMA solution to solar anomaly is shown in 
fig. 3. 

Such small neutrino degeneracy cannot have a noticeable cosmological impact, in 
particular, on LSS and CMBR. However, bound ( 1391) can be relaxed [TH], if neutrinos 
have a new interaction with light scalar bosons, Majorons. The potential of neutrinos 
induced by the Majoron exchange inhibits early oscillations in the plasma. 




Figure 3. Transformation of initial chemical potential of into those of v e and v T as a 
function of the plasma temperature, with and without neutrino self-interactions. 



8 Cosmological impact of neutrino statistics viola- 
tion 

According to the famous Pauli theorem, particles with integer spins obey Bose-Einstein 
statistics, while those with half integer spin obey Fermi-Dirac statistics. There is no con- 
sistent formulation of a theory where this theorem is violated. Nevertheless phenomeno- 
logical manifestations of spin-statistics violation are discussed and experimental bounds 
are presented. Such bounds are very strong for electrons and nucleons, while for neutrinos 
a large spin-statistics violation is allowed. Such a possibility is discussed in ref. [20]. In 
this paper one can also find a list of references to theoretical and experimental works on 
the violation of the Pauli theorem. 

The best experimental bound on breaking of the exclusion principle for neutrinos 
follows from the double beta decay but still even this best bound allows for about 50% 
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violation [21]. In other words, in principle neutrino may be half and half of fermionic and 
bosonic parts. In fact, the existing data allow neutrinos to be more than half bosonic, up 
to 60-70%. 

If neutrino statistics is not purely the Dirac-Fermi one, there may be very interesting 
cosmological consequences. Evidently the theoretical predictions for abundances of light 
elements created at BBN would change [22J. If neutrinos were purely bosonic, their energy 
density would be 8/7 of normal fermionic u, giving AN U = 3/7. On the other hand, a 
larger number density of v e would lead to smaller temperature of n/p-freezing and thus 
is equivalent to negative AN U . The net result is that the effective number of neutrinos at 
BBN should be smaller than three: Njff^ = 2.43. This is exactly the number obtained 
from the data analysis in ref. |18j . 

In ref. [22] it was argued that the equilibrium neutrino distribution for mixed statistics 
has the form: 

= [ex V (E/T) + Kr 1 . (40) 

where the parameter k changes from (-1) to (+1) and interpolates between Bose and 
Fermi statistics. 

Another interesting effect of (partly) bosonic neutrinos is that they might form cos- 
mological Bose condensate. In the (excluded) case of purely bosonic neutrinos their 
equilibrium distribution would be 

U = U~P 1 WT 11 + OS^(k), ( 41 ) 

exp[(E - fi v )/T - 1\ 

where the chemical potential reaches the maximum allowed value for bosons: /^ max ) = m u . 
If the equilibrium distribution of neutrinos is mixed one (1401) . with (— 1) < k < the 
maximum value of the chemical potential, at which Bose condensation of neutrinos could 
take place, would be y^ aax ^ > = m v — Tin \ k\. 

If neutrinos might Bose condense then they can make not only hot dark matter but 
also (all) cold dark matter and we do not need for that any new particles or other objects. 
One can avoid the lower bound on neutrino mass [23J based either on Fermi statistics or 
(for bosons) on the Liuville theorem [21] because of the condensate formation. So the 
cosmological dark matter can consist of old known particles (massive, though very light 
neutrinos) but to this end we need very new physics. 

If the neutrino condensate makes all cosmological dark matter, its average number 
density should be about 10 4 (m I/ /0.1 eV)/cm 3 , while in the Galaxy it can be 5-6 orders of 
magnitude larger. Probably it is not enough to explain the observation of electrons out 
of the kinematically allowed region in the tritium beta decay, mentioned at the end of 
sec. El 



9 Neutrino oscillations in the early universe 

Here we will consider possible effects on BBN if the usual active neutrinos are mixed with 
new sterile ones. Neutrino propagation in plasma is governed by the effective potential 
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(refraction index) [25] : 

V e a ff = ±C lV G F T 3 + cf^^- , (42) 

where the constants Cj are of the order unity and 77 is the plasma charge asymmetry: 

r/ £ ) = 2rj Ue + r) Vfl + t] Ut +rj e - rj n /2 (for v e ) (43) 

V M = 2r]^ + r] Ve + r) VT - r] n /2 (for (44) 

The first term comes from the averaging of the neutrino current ( J ) over plasma. It has 
different signs for neutrinos and antineutrinos. It is proportional to the charge asymme- 
try in plasma, i.e. to the difference of the number densities of particles and antiparticles 
and in cosmological situation it is usually small but in stars or in the Earth it gives the 
dominant contribution. The second term originates from the non-locality of weak inter- 
actions related to the exchange of the intermediate bosons. Because of that neutrino can 
be absorbed in one space-time point and re-emitted in different one. This term is usually 
negligible in astrophysics but is dominant in cosmology especially at high temperatures, 
T ~ 10 MeV. 

Mixing between active and sterile neutrinos would lead to the following phenomena 
which in turn could produce noticeable effects on BBN: 

1. An increase of the number of particle species due to production of sterile neutrinos. It 
leads to > 3. 

2. Since the probability of the oscillations depends upon neutrino energy, the neutrino 
spectrum could be different from the usual equilibrium one. As we mentioned above, BBN 
is very sensitive to the spectrum distortion of v e . 

3. Generation of large lepton asymmetry in the active neutrino and especially in v e sector 
by the resonance oscillations between active and sterile neutrinos [26] . 

For more details see review [1] 

In the simplified case when the mixings of active neutrinos are neglected, and the 
mass difference between active and sterile neutrinos is positive, so there is no resonance 
transition, the bounds on the mixing angle and the mass difference between one active 
and one sterile neutrino are: 

(5ml eU JeV 2 ) sin 4 28 UeU ° = 3.16 • 10" 5 ln 2 (l - AN V ) (45) 
{5ml^JeV 2 ) sin 4 26^ Us = 1.74 • 1(T 5 ln 2 (l - AN U ) (46) 

The impact of oscillations between active and sterile neutrinos for realistic mixing 
between active neutrinos for both resonance and non-resonance cases was considered in 
ref. [2?]. Unfortunately the results cannot be presented in a simple analytical form, as 
above, but only as figures. 

10 BBN and right-handed neutrinos, Wr, Zr, neu- 
trino magnetic moments, etc 

We mentioned above that if neutrino mass is non-zero, there must be both helicity states, 
right-handed and left-handed. Now let us estimate the rate of production of the "wrong" 
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right-handed neutrinos in the early universe. The most favorable period for the production 
of vr in the standard theory, i.e. without right handed currents, only by neutrino mass, 
took place at the temperatures of the order of W and/or Z mass, when these intermediate 
bosons were present in the plasma. The production rate of vr by the decays of W and Z 
is equal to: 

Tr = ^;~ 10 {~t) — — ' (47) 

where T\y Z is the decay width of W or Z into neutrino channel and nw,z are the number 
densities of these bosons. Comparing with the expansion rate, H ~ T 2 /mpi, one sees 
that vr were never produced abundantly if m u respected the GZ-bound. 

If in addition to the usual weak intermediate bosons, Wl and Zl (we omitted above the 
sub-index L) there exist right-handed bosons which interact with right-handed fermions, 
the probability of interactions of e.g. right-handed neutrinos would be inversely propor- 
tional to m^ R . The energy density of v R at BBN should be smaller than, say, 0.3 of the 
energy density of the usual left-handed neutrinos. It demands that u R should decouple 
before the QCD phase transition. If so, the right-handed neutrinos would be diluted to 
the necessary amount by the entropy released in the massive particles annihilation. This 
condition leads to the lower bound of W R mass: 

^>2.5TeV (^^-) 4/3 (48) 
m WL V200MeV/ v ; 

This is an order of magnitude better than the direct experimental limit. 

There is a different model for generation of right-handed neutrinos by the mixing 
between W R and Wl- 

W x = cos 9W L + sin 9W R (49) 
In this case the production of vr at T < Tqcd is given by 

r = r R /tf = sin 2 4%^ 3 (50) 



2MeV7 

Demanding r < 0.3 results in sin 2 9 < 10~ 6 . It is interesting that the effect does not 
vanish for — > oo. 

The increased total number of neutrino species, which would lead to a larger value of 
(n/p)-ratio, may be compensated by an increase of the number of the right-handed elec- 
tronic neutrinos which makes the n — p transformation more efficient and thus diminish 
their frozen ratio. However the interaction of vr with nucleons is proportional to ~ sin 2 6 
and the effect is small, if 9 < 10 -3 , as is written above. On the other hand, a large 9 ~ 1 
may be allowed because in this case the compensation would be possible. 

If neutrinos have a non-zero magnetic moment, \i v , it also leads to production of 
right-handed helicity states. For example, vr could be produced in the electromagnetic 
reaction: 

e ± + v L -> e ± + v R (51) 
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Demanding that right-handed neutrinos contributes not more than AN U = 0.5, we obtain: 



(52) 



where fig is the Bohr magneton. 

If there exists primordial magnetic field, coherent on sufficiently large distances, right- 
handed neutrinos can be produced by spin flip in this field. In this case the magnetic 
moment is bounded by: 



The magnitude of the B primor d is unknown but if we assume that the intergalactic magnetic 
fields have the strength B int _ ga i ~ 10 -6 Gauss and that they were generated in the early 
universe and evolved adiabatically to the present time, their strength at T = 1 MeV 
should be 10 13 G and we obtain a very strong bound: 



11 Conclusion 

Thus we see that cosmology and astronomy happen to be the most efficient neutrino 
"detector" at the present time. The upper bound on neutrino mass at the level of a fraction 
of eV is noticeably stronger than that obtained by direct experiments. On the other hand, 
the cosmological bound on m v may be somewhat relaxed if new interactions of neutrinos 
exist, which led to their enhanced annihilation in the early universe and diminished the 
neutrino number density, or there is a new light long-lived particle whose electromagnetic 
decay might dilute the neutrino-to-photon ratio. Moreover, a modification of the spectrum 
of the cosmological density perturbations at a few Mpc scale would invalidate the most 
restrictive cosmological bounds on m v . So KATRIN would be helpful to diminish the 
level of ambiguity. 

The number of neutrino species is now best restricted by BBN, N v = 3 ± 0.3. The 
angular spectrum of CMBR is less efficient, giving N v — 3 ± 1. But one can hope that 
the Planck mission will be competitive with BBN or even will be better than that. One 
should keep in mind, however, that the number of the effective neutrino species found 
from BBN and CMBR are sensitive to physics at different cosmological periods and they 
are not not necessarily to coincide. 

The heating of neutrinos by hotter electrons and positrons, when the universe was 
about 1 sec old, distorted the spectrum of neutrinos and increased their energy density 
by 3%, plus 1% from plasma corrections. If this effect is observed in the future it will be 
the direct measurement of a physical process which took place in the very early universe. 

The bound on the cosmological lepton asymmetry for all neutrino flavors found from 
BBN, was significantly improved when it was established that the active neutrinos are 
strongly mixed. The maximum allowed value for the leptonic chemical potential, fi/T < 
0.07, makes it cosmologically insignificant. 

The bounds on the mixing between active and possible sterile neutrinos are signifi- 
cantly modified because of the strong mixing of all active neutrinos between themselves 




(53) 



Hu < 10 a Vb 



(54) 
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in comparison with earlier calculations, when the mixing of active ones were not taken 
into account. Sometimes the new bounds are more restrictive than the old approximate 
ones and in all cases are better than direct experiment results. 

Both BBN and LSS formation are sensitive to violation of Fermi statistics for neutrinos. 
If the bosonic fraction in neutrinos is larger than 50%, cosmic neutrinos may condence 
and make all dark matter in the universe. Unfortunately there is no consistent theory of 
violation of the spin-statistics theorem and the consequences of its violation are not well 
justified. 

Right handed neutrinos which must exist, if m v 7^ 0, are not allowed to be abundant 
at BBN and one can derive from that quite strong bounds on the mass of possible right- 
handed intermediate bosons and mixing of them with the left-handed ones. Similarly the 
bound on magnetic moment of neutrinos can be derived. All these bounds are stronger 
than those found from direct experiment. 

Acknowledgement. I thank Gary Steigman for critical comments. 

References 

[1] A.D. Dolgov, Phys. Rept. 370 (2002) 333. 

[2] S.S. Gerstein, Ya.B. Zeldovich, Pis'ma Zh. Eksp. Teor. Fiz. 4 (1966) 174; English 
translation: JETP Lett. 4 (1966) 120. 

[3] R. Cowsik, J. McClelland, Phys. Rev. Lett. 29 (1972) 669. 

[4] W.-M. Yao et al., Journal of Physics, G 33 (2006) 1. 

[5] G.F. Giudice, E.W. Kolb, A. Riotto, Phys. Rev. D64 (2001) 023508. 

[6] Y. Farzan, Phys. Rev. D67 (2003) 073015. 

[7] A.D. Dolgov, M. Fukugita, Pisma Zh. Eksp. Teor. Fiz. 56 (1992) 129; [JETP Lett. 56 
(1992) 123]; 

A.D. Dolgov, M. Fukugita, Phys. Rev. D46 (1992) 5378; 
S. Dodelson, M.S. Turner, Phys. Rev. D46 (1992) 3372. 

[8] R.E. Lopez, S. Dodelson, A. Heckler, M.S. Turner, Phys. Rev. Lett.82 (1999) 3952. 

[9] K. Ichikawa, M. Fukugita, M. Kawasaki, Phys. Rev. D71 (2005) 043001; 

M. Fukugita, K. Ichikawa, M. Kawasaki, O. Lahav, Phys. Rev. D74 (2006) 027302. 

[10] J. Lesgourgues, S. Pastor, Phys. Rept. 429 (2006) 307; 
S. Hannestad, Ann. Rev. Nucl. Part. Sci 56 (2006) 137. 

[11] S. Hannestad, G.G. Raffelt, JCAP 0611 (2006) 016. 

[12] D. Fargion, B. Mele, A. Salis, Astrophys. J. 517 (1999) 725; 
T.J. Weiler, Astropart. Phys. 11 (1999) 303. 



20 



[13] V.F.Shvartsman, V.B.Braginskii, S.S.Gershtein, Ya.B.Zeldovich, M.Yu.Khlopov, 
Pis'ma v ZhETF 36 (1982) 224 [English translation: JETP Lett. 36 (1982) 277]. 



[14 
[15 

[16 
[17 
[18 
[19 
[20 
[21 

[22 
[23 

[24; 

[25 
[26 



S. Weinberg, Phys. Rev. 128 (1962) 1457. 

A.D. Dolgov, S.H. Hansen, S. Pastor, S.T. Petcov, G.G. Raffelt, D.V. Semikoz, Nucl. 
Phys. B632 (2002) 363. 

A.G. Cocco, G. Mangano, M. Messina, larXiv:0711.1762l [astro-ph]. 

V.M. Lobashev, Nucl. Phys. Proc. Suppl. 91 (2001) 280. 

V. Simha, G. Steigman, ilarXiv: 0803 .34651 

A.D. Dolgov, F. Takahashi, Nucl. Phys. B688 (2004) 189. 



A.D. Dolgov, A.Yu. Smirnov, Phys. Lett. B621 (2005) 1; arXiv:hep-ph/0501066 



A.S. Barabash, A.D. Dolgov, R. Dvornicky, F. Simkovic, A.Yu. Smirnov, Nucl. Phys. 
B783 (2007) 90: ! arXiv:0704.2944l [hep-ph]. 

L. Cucurull, J. A. Grifols, R. Toldra, Astropart. Phys. 4 (1996) 391, |astro-ph/9506040[ 
A. D. Dolgov, S. H. Hansen, A. Yu. Smirnov, JCAP 0506 (2005) 004, 
|astro-ph/05036T2l 

S. Tremaine, J.E. Gunn, Phys. Rev. Lett. 42 (1979) 407. 

J. Madsen, Phys. Rev. Lett. 64 (1990) 2744; Phys. Rev. D 44 (1991) 999. 

D. Notzold, G. Raffelt, Nucl. Phys. B307 (1988) 924. 

R. Foot, M.J. Thomson, R.R. Volkas, Phys. Rev. D53 (1996) 5349; 
P. Di Bari, R. Foot, Phys. Rev. D65 (2002) 045003. 
A.D. Dolgov, Nucl. Phys. B610 (2001) 411. 



[27] A.D. Dolgov, F.L. Villante, Nucl. Phys. B679 (2004) 261. 



21 



